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A systematic theory of the conductance measurements of non-invasive scanning gate microscopy is 
presented that provides an interpretation of what precisely is being measured. A scattering approach 
is used to derive exphcit expressions for the first and second order conductance changes due to the 
perturbation by the tip potential in terms of the scattering states of the unperturbed structure. In 
the case of a quantum point contact, the first order correction dominates at the conductance steps 
and vanishes on the plateaus where the second order term dominates. Both corrections are non-local 
for a generic structure. Only in special cases, such as that of a centrally symmetric quantum point 
contact in the conductance quantization regime, can the second order correction be unambiguously 
related with the local current density. In the case of an abrupt quantum point contact we are 
able to obtain analytic expressions for the scattering eigenfunctions and thus evaluate the resulting 
conductance corrections. 



I. INTRODUCTION 

Understanding the underlying mechanisms of the elec- 
tronic transport properties of nanostructures in the quan- 
tum coherent regime is of great fundamental interest, 
and also of foremost importance for possible applica- 
tions in nanoelectronics devices and quantum comput- 
ing. The scanning gate microscopy (SGM) technique, 
born more than a decade ago,'^' provides additional infor- 
mation about coherent transport beyond that obtained 
in traditional transport experiments, and has therefore 
attracted considerable interest. 

In SGM, the transport through a nanostructure is mea- 
sured while the charged tip of an atomic force microscope 
(AFM) is scanned over the sample. The tip scatters elec- 
trons and thereby influences the sample's transport prop- 
erties. One then studies how transport coefficients, for 
example the linear conductance through the nanostruc- 
ture, depend on the position of the tip over the sam- 
ple. The technique (for a brief review see Ref. has 
been applied to quantum point contacts (QPCs) defined 
in two-dimensional electron gases (2DEGs),'2Kil and to a 
variety of other systems, including carbon nanotubeSj^i^l 
quantum dots fabri cated in various systems (semicon- 
ductor 2DEGs,ll^'J^ carbon nanotubes,!^ semiconductor 
nanowires,^^ and graphene!^, small Aharon ov-Bohm 
rings, ""^^ ^'^ Hall bars,^^ and edge channels .^221 

Recently, SGM has been extended to so-called scan- 
ning probe microscopj^, where the scanning of a 
charged tip together with the measurements of the re- 
sulting Coulomb blockade peaks energy-shifts of a one- 
dimensional quantum dot enables the extraction of the 
electronic wave- functions density profile. 

Based on the plausible argument that the AFM tip 
has the strongest infiuence on the conductance if it is lo- 
cated directly above a region of large current density, it is 
tempting to interpret these conductance maps in terms 
of the local "electron fiow". This was proposed in the 



first SGM measurements on quantum point contacts,'^ 
where branching structures were observed. For QPC- 
tip distances below the phase coherence length, interfer- 
ence patterns were observed as well. Since they extended 
to regions more than a thermal length away from the 
QPC, those patterns were attributed to the interference 
between electron paths scattered back to the QPC by 
the tip, and by an impurity within a thermal length of 
the tip.l^il Such interference patterns disappear for very 
clean sampled except at very low temperatures,^ where 
the thermal length exceeds the QPC-tip distance. In- 
terestingly, the pattern for a QPC tuned to the first con- 
ductance plateau has been observed to be rather different 
from the one obtained when the QPC is tuned to the first 
step.f^ The spacial periodicity of the interferences when 
the tip moves away from the QPC was expected to be 
half a Fermi wave-length, assuming it would be due to 
the interference between the electron path scattered back 
by the tip with another one that does not visit the tip. 
Though this is consistent with most experimental obser- 
vations, recent experiments show different periodicity in 
different regions of the same sampld^ and challenge this 
simple picture. 

While the interpretation of SGM data in terms of elec- 
tron flow seems consistent with numerical studies of the 
local current density and the SGM response in QPCd^SESlj 
SGM experiments on small ring structures combined with 
numerical simulation j^^l^^l have concluded that in this sit- 
uation the conductance change is related to the local den- 
sity of states rather than to the local current. Though 
a generally applicable theory leading to an unambigu- 
ous and quantitative interpretation of the SGM data is 
still lacking, experiments continue to propose novel and 
interesting uses of the SGM technique. For example, 
while theoretical modeling suggested looking for inter- 
action effects in the regime where the tip is very close to 
the QPC, experimental evidence indicates that SGM 
in the nonlinear transport voltage regime gives access to 
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information about electron-electron scattering inside the 
sample.® 

In spite of the numerous measurements to be found in 
the literature, it is not always obvious how the now con- 
siderable quantity of data is to be interpreted. Further 
theoretical studies of coherent transport in the presence 
of a local scatterer are needed to fill in the gaps of current 
understanding. 

In a previous paper, we initiated a systematic ap- 
proach to the theory of the conductance change induced 
by a local perturbation. The resulting general expres- 
sions allow one to calculate the correction to the con- 
ductance to first and second order in the perturbation 
caused by the tip, starting from the scattering properties 
and wave-functions of the unperturbed nanostructure. It 
was found that the first order correction is suppressed 
when the QPC is tuned to a conductance plateau. This 
explains a change in the SGM response depending on 
the tuning of the QPC, and is consistent with a thermal 
enhancement of the interference patterns proposed for a 
QPC tuned close to the edge of a plateaii^. 

This paper provides a detailed derivation of the expres- 
sions presented in Ref . pS] , and addresses the question of 
the relationship between the SGM signal and local quan- 
tities, such as the current density, by applying the general 
formulas to the case of the SGM response of a QPC with 
abrupt openings. This model is capable of reproducing 
quite accurately and in a controlled way the features re- 
lated to the QPC conductance quantization, and, cru- 
cially, allows one to obtain expressions for the scattering 
wave-functions, thus making it possible to compare the 
SGM response to the local current density. These two 
quantities are directly related only if: (i) the QPC is 
tuned to a conductance plateau, and (ii) the symmetry 
of the structure is such that the wave leaving the QPC 
can be written as the product of a radial and an angu- 
lar part, as it is for an outgoing wave from a point-like 
source in a clean medium. However, when the symmetry 
of the structure is broken by imperfections of the sam- 
ple or by disorder in the neighboring 2DEG, even on a 
plateau the SGM response is no longer given by the local 
current density. 

The paper is organized as follows. Section |TT] presents 
the scattering formalism used to evaluate the SGM re- 
sponse. The derivation of the first and second order con- 
ductance corrections follows in Sees. |III| and |IV| respec- 
tively. The calculation of the scattering wave-functions 
and the scattering eigenmodes for an abrupt QPC is pre- 
sented in Sees. |V] and |VI[ The symmetry of the unper- 
turbed structure is shown in Sec. |VII| to play a key role 
for the determination of the conductance corrections, and 
the example of an abrupt QPC is used for the discussion 
of the connection between conductance corrections and 
current densities. A summary is presented in the con- 



cluding Sec. VIII Details of the energy integrations used 
in Sees. |IIl| and |IV| alternative routes to reach the main 
results, the comparison of our expressions with exact re- 
sults available for one-dimensional models, and that with 



the formalism proposed in Ref. 
pendices. 



all appear in the ap- 



II. MODEL DESCRIPTION AND SCATTERING 
FORMALISM 

The conduction electrons in a typical SGM setup are 
described by a Hamiltonian 



H = Hn 



Vt 



(1) 



where Hq represents the unperturbed structure to be 
characterized and Vt the electrostatic potential gener- 
ated by the perturbing tip. While most of the analysis 
is developed for an arbitrary Hq , the case of a quantum 
point contact is treated in detail. Similarly, general re- 
sults for a non-specified Vt are established, and then the 
case of local spatial perturbations is the focus. 

The scattering theory of quantum conductanc^^^^'Sll as- 
sumes that the unperturbed structure, considered as a 
non-interacting scatterer, is connected to electron reser- 
voirs through disorder-free leads of finite cross section 
and semi-infinite in the longitudinal direction. Even 
though the scattering theory is applicable to an arbitrary 
number of leads and spatial dimensions, we will restrict 
ourselves to the case of interest of a two-lead setup within 
a two-dimensional space spanned by vectors r = {x,y). 

The electrons in the reservoirs are assumed to be free, 
and therefore their dispersion relation is e = h'^k'^ /2Mc, 
with k the magnitude of the two-dimensional wave- vector 
and Afc the effective electron mass. Taking the x- 
direction as the longitudinal one, the incoming (— ) and 
outgoing (-1-) modes in lead 1 (left) and 2 (right) with 
energy e can be taken, respectively, as 



exp[±ifcja;] </)a(y) , a; < , (2a) 



exp[=FifcJa;] (j)a{v) , x> Q . (2b) 



With 0a (y) we denote the wave- function of the a trans- 
verse channel (with quantized energy ei'^), whereas ka 
is the longitudinal wave-vector, Va — hka/M^ the longi- 
tudinal velocity, and Ea^ ~ fi^fc^/2Afc the longitudinal 
energy. The relationship determines, for 

each lead, the N propagating modes which satisfy fc^ > 
at energy e. We always take ka > 0. The notation k^ 
stands for an infinitesimal negative (positive) imaginary 
part given to ka for incoming (outgoing) modes. Choos- 
ing the constant c equal to 1 amounts to the so-called 
unit-fiux normalization for the lead modes. We adopt 
however a slightly different convention and set 



(3) 



This is notationally simpler and leads to a current den- 
sity in the x-direction and a current (per spin and unit 
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energy) associated with the right-(left-)moving mode 
1(2), e, a, respectively given by ±e / h\(j)a{y)\'^ and ±e/h. 
The overall signs result from our convention of taking as 
positive the current of positive charges moving from left 
to right. 

Although any separable potential can be treated if al- 
lowance is made for an a;-dependence of ka , we consider a 
confining potential that is x-independent in the asymp- 
totic regions. For simplicity we choose a hard wall con- 
finement in the y direction by taking leads of width 2W. 
Then 



.(t) _ 



2Mn 



and 



My) 



(-IF 



sin[(ja(y- VK)] , 



(4) 



(5) 



where qa = TTa/2W is the transverse wave- vector satisfy- 
ing ka — y^k"^ — Qai and p = Int{a/2}. The introduction 
of the phase (—1)'' is a matter of convenience for writ- 
ing the scattering wave-functions, as done in Sec. |IV[ It 
merely gives a sign alternation within each of the families 
of even and odd (in y) modes corresponding, respectively, 
to odd and even a. 



Once a quantum coherent scatterer (of linear extension 
L in the x direction) is placed at the coordinate origin, 

the incoming modes ^^2) e a S^^^ ^'^^^ outgoing scat- 
tering states (defined for all x) that in the asymptotic 
regions are, respectively. 




(-) 

1,6, a 
N 



Etina^t\i^), x«-L/2 



X > L/2 



(+) 



(+) 



X < -L/2 



(6a) 
(6b) 



In order to simplify the notation we do not write the stan- 
dard + label corresponding to outgoing scattering states. 
(The incoming scattering states generated from the out- 
going modes ^^2) e a '^^^ used in this work.) The 
index (0) is used for labeling unperturbed quantities de- 
pending on Hq only. The N x N matrices r (r') and t (t') 
characterize, respectively, the reflection and transmission 
matrices from lead 1 (2). The 2Nx2N scattering matrix 
S, relating incoming and outgoing modes, is given by 



r t' 
t r' 



(7) 



Current conservation dictates the unitarity of the scat- 
tering matrix {SS^ — I)- In the absence of magnetic 
fields, which is the situation considered in this work. 



time-reversal symmetry implies that 5* is a symmetric 
matrix {S^ — S). For simplicity we do not explicitly 
write the energy dependence of the various components 
of the scattering matrix, nor do we include the index (0) 
that should in principle be assigned for consistency of the 
notation. 

The normalization (3]) ensures that the scattering 
states constitute an orthonormal basis verifying 



dr ^\%{r) m^ir) ^ Sa Sis - e) Sa 



,(0) 



(8) 



and have the same flux normalization as the lead states. 

The transmission and reflection amplitudes between 
modes a and b can be obtained, respectively, from the 
retarded Green function C/'^"^ (r, f , e) associated to the un- 
perturbed Hamiltonian Hq [3T] 



t,a - iHvaVi,y/^ exp [-t{k+x-k+x)] / dy dy 4>l{y) g^''\v,v,e) Mv) , M 
rba = -Sab exp [i{k+x + k^x)] exp - x\] 

+ ihivaVby^^ e^p[-t{k+x + k+x)] f dy f dy ^.iy) g^°Hr,v,e) Mv) ■ (9b) 



The integrations take place at the transverse cross sec- tions Sx on the left lead and Sx on the right (left) lead for 
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the transmission (reflection) amplitudes. Since only re- 
tarded Green functions are used, the standard (+) label 
for Q^^^ is dropped. 

The Landauer-Biittiker approach to quantum trans- 
port derives the conductance of a coherent scatterer tak- 
ing as building blocks the current carried by the scat- 
tering states,- ^ in accordance with the fact the in the 
traditional conductance measurement setups only the to- 
tal integrated currents are relevant. However, the SGM 
technique probes the sensitivity of the current to local 
perturbations, thus providing a spatial resolution that 



yields information about the scatterer beyond that of the 
traditional setup. In order to understand the outcome of 
SGM measurements, and the perturbative approach that 
is developed, it is useful to reformulate the scattering 
formalism starting from the current operator, defined as 



2Mr 



[p5(f-r) + (5(f-r)p], 



(10) 



where f and p denote the position and momentum op- 
erators, respectively. The matrix elements of the x- 
component of the current density in the scattering states 
basis read 



eh 
2iAL 



l.£,a 



dx ''^^5 



(r) 



(11) 



The diagonal matrix element, that we write j^i^f^i^), rep- 
resents the current density per spin and unit energy as- 
sociated with the state ^f'} g^- For a given incoming lead 
I and energy e it is useful to define an N x N current ma- 
trix I;''^' whose elements in the scattering states subspace 
I, e are given by 



AO) 



dy r(°)(r) _ (e,£) 



(12) 



The independence of from the cross section Sx cho- 
sen for the integration is a consequence of current con- 
servation. Given the one-to-one correspondence between 
incoming modes and outgoing scattering states, and the 
asymptotic form (|6| of the latter, the current matrix ele- 
ments (involving scattering states) can be identified with 
those of tH (involving lead modes), that is. 



(0) 



(13) 



The diagonal matrix element T-^x\ is the current (per 

L ' J aa 

spin and unit energy) associated with the scattering state 
1, e, a 



N / ^ \ 

6=1 \ 6=1 / 



(14) 



The total current from left to right can be written as 



AO) _ 
^1 — 



Ml 



N 



(15) 



with Pais) — (j^hva)^^ the one-dimensional density of 
lead modes (including the spin degeneracy factor) and 
P-i(2) the electrochemical potential in the left (right) 



Working in linear response to the applied bias eV = 
fii—fi2 leads to the two-probe Landauer-Biittiker formula 
for the dimensionless conductance of the scattereil^ 
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(0) 



1 / 



(0) 



2e^/h V 



(16) 



where all quantities are calculated at the Fermi energy ep 
of the reservoirs, and the trace is taken over the incoming, 
right-moving modes. 

Within the approach, it will be convenient to 
write the sy mmetr ic scattering matrix S in the polar 
decompositioiP^I^ 



S = 



T 









-n 
r 



T 

n 



ui 

U2 



(17) 



where Ux(2) are N x N unitary matrices, while TZ and 
T are diagonal matrices defined by the reflection and 
transmission eigenvalues, respectively given by 



I rn 



1 

1 + A, 



1/2 



1/2 



(18a) 
(18b) 



with Am real positive. 

Since tH = u\T^ui and t'H' — 1/2X^1*2, the matrices 
Mi(2) diagonalize, respectively, tH and t'H'. The trans- 
mission eigenvectors, or transmission eigenmodes, are of 
the form 



N 



9i 



(-) 



(r)=E["i]™a ^UaW , a;< 0, (19a) 



a=l 
N 



9i:J,n.ir) = E ^t]air) , X > . (19b) 



reservoir. 
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The dimensionless conductance only depends of the 
transmission eigenvalues Tm as 



N 

,(0) ^ ^ 



(20) 



Our goal is to determine the change of the conductance 
from (7*^°^ to g when the full Hamiltonian H of Eq. ([I]) 
is considered. In the perturbative approach in powers 
of Vt developed in Sees. |III| and |IV| the properties of 
the unperturbed system are important, particularly the 
eigenvalues and eigenfunctions of the current matrix x[^^ . 
In Sees. [V| and |VI| we will focus on QPCs. building such 



eigenfunctions and then studying in Sec. VII the relation 
between the unperturbed current densities and the tip- 
induced conductance changes. 



III. FIRST-ORDER PERTURBATION IN THE 
TIP POTENTIAL 



In order to analyze the effect of the tip voltage in 
its least invasive form we will treat the potential Vt 
of Eq. as a perturbation to Hq. Two substantially 
equivalent approaches are possible. One can compute 
the corrections to the retarded Green function of the un- 
perturbed structure G^^^ via the Dyson equation, thus 
obtaining the modified transmission and refiection ampli- 
tudes through Eqs. ([9|. Alternatively, the corrections to 

the unperturbed scattering states ^'|*'j'a(r) can be com- 
puted by means of the Lippmann-Schwinger equation. 



and the modified states can then be employed to obtain 
the corresponding current density and full current. We 
mentioned that the computation of the local current den- 
sities is fundamental for our purposes, as these need be 
compared to the SGM-induced conductance corrections 
in order to establish the relation - if any - between the 
two quantities. We will therefore follow the second route, 
pursuing further the reasoning of the previous Section. A 
discussion of the Green function approach is given how- 
ever in Appendix |D] The Lippmann-Schwinger equation 
relates the unperturbed scattering state ^'l^^'' ^ (r) to the 
corresponding perturbed state ^'/,e^a(r) according to 

*^.,a(r) - *W Jr) +Jdr ^(r, r, £)^T(f)*/,e.a(f) . 

(21) 

The Green function Q^^\ introduced in Eq. has the 
following spectral decomposition in the scattering states 
basis 



2 .oo ^- N 



-e , 

1=1 a=l 



(22) 

having defined = e ± ir], with r] positive infinitesimal. 

In Eq. (21 1 the first order correction in Vt (Born ap- 
proximation) to the scattering state is obtained as 

*Sa(r) = /df g(")(r,f,£)l/T(f)<)„(f) . (23) 

The corresponding change of the current density associ- 
ated with the state l,e,a is given by 



€!aW-25:Re 



N 



(t) £ 



+ 



1=1 1 ""i s=i 
where the matrix elements of the perturbing potential in the scattering states basis are 



(24) 



(25) 



The change of the current associated with the scattering state l,e,a, obtained by integrating Eq. (24 1 over a cross 
section Sx in the right lead, is 



eh 
M 



N 



(26) 



The final result is obviously independent of x and the chosen lead (by current conservation) . In Eq. ( 26 1 we have 
introduced a shorthand notation for quantities involving the unperturbed states 



JV 



6=1 



(27a) 



N 



^ ' ^L^baexp [i(fc+ - fcj )x\ \ , 



(27b) 
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where the barred longitudinal wave-vectors and scat- 
tering amplitudes are taken at the total energy e, and 
N — minjA'', N}. Performing the e integration by con- 
tour in the complex plane (see Appendix |A| leads to 



with the matrix I^^j given by 



(1) 



-^1 



(1) 



(28) 



= ^Imji^t V" (£,£)+ (£,£)} . (29) 



The Nx N matrix V" (£,£), spanning the space of incom- 
ing modes 1,1, has its element (a, a) defined by Eq. (25 1 



This definition is based in the one-to-one correspondence 
between incoming modes and outgoing scattering states. 
Using the unitarity of S and the fact that V^''(j(£,£) is a 
real quantity we can write Eq. ( 29 I as 



= -^lia{r^r V^\e,e) +rh' V^\e,e)} . (30) 

The change of the total current, up to first order in Vi- 
and in linear response to the applied voltage V, is ob- 
tained [as in Eq. ( 15 )] by summing the contribution of 
all modes at the Fermi energy. The first term of (30) 



vanishes when summed over a and therefore the 0{Vt 
change of the zero-temperature dimensionless conduc 
tance i^ 



-4t7 Im {Tr [rh' V^^] } , 



(31) 



where all quantities are evaluated at £p and the trace 
is over the incoming, right-moving modes. This result 
is valid for the general situation where quantum trans- 
port through a scatterer is modified by a weak perturba- 
tion. The matrix rH' depends only on the unperturbed 
scatterer, while the tip's effect appears in the V matrix 
elements. 

The conductance change is invariant under a change of 
current direction as well as under the changes Vt{x, y) — 
Vt(±x, ±y) of the perturbing potential, provided the un- 
perturbed system verifies the corresponding (up-down 
and/or left-right) refection symmetries. 



The conductance's sensitivity to electrostatic potential 
variations was considered in a one-dimensional geometry 
in Ref . 1301 The example of a (5-function barrier perturbed 
by a local tip can be analytically calculated. We show 
in Appendix [B] that this particular case is in agreement 
with the general expression (31 1. 

Like g^'^\ the conductance correction g'^^^ is given by 
the trace over the space of propagating modes in one 
lead. The obvious requirement of a base-independent re- 
sult is thus fulfilled. We stress that, even if the matrix 
V^^ is obtained from the overlaps of scattering states, 
the entries of its matrix elements are mode indices. The 
relationship between scattering states and lead modes is 
fixed by Eq. ([6|, and a base transformation in the space 
of the modes of lead 1 (2) induces a change of '^i,e,a and 
^2,e,a for X < (a; > 0). 



Using the polar decomposition (17) of the scattering 
matrix we can write 
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(1) = 47r 



N 



TZraTra ItH [UH^] 



(32) 



Where Z^"(£, e) = U2V"(£, s)u\ represents the matrix 
upon transformation into the transmission eigenmodes 
basis (19). In Eq. (32), as in Eq. (31), the energy argu- 



ments are understood to be taken at £f. 

From Eq. ( 32 I we see that the m'^ mode contribution 
to the conductance correction g'^^ is appreciable only if 
m is partially opened. This observation is of foremost 
importance in the case of QPCs, where the transmission 
eigenmodes open one by one. On the w}^ conductance 
plateau we have Tm — 1 and TZm — 0, therefore g^^-* 
is suppressed. The first order conductance correction is 
only relevant in the vicinity of a conductance step, and 
only the contribution arising from the transmission eigen- 
mode that is partially open matters. In order to capture 
the dominant correction on the conductance plateaus, it 
is necessary to go beyond the first order Born approxi- 
mation. 



IV. SECOND-ORDER PERTURBATION IN THE 
TIP POTENTIAL 

The second order correction in Vr to the scattering 
state of Eq. (21 1 is obtained as 



(33) 



The second-order change of the current density associated with the scattering state l,e,a is 



(34) 
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with 



eh 
Ml 



Im|vl/(o)*(r)Ax],p) M-vI/p) (r)Avl>W* (r) 

I l,e,a\ ) l,e,a\ I l,e,a\ ) l,e,a\'- ) 



2^ Re 

1,1=1 
eh 

2 



d£ 



de 



AT N 



(') e+ — £ /.(') e+ — £ 



El-' 



1.1=1 



d£ 



(t) £' 



a=l a=l 



W AT 



(35a) 



£ ./^(t) £+ — £ 



EE[^]i(=>^)*S;w|*g,(r)[^]L(^,e)L (35b) 



The calculation of the corresponding current corrections follows the lines presented in the last section, but it is 
considerably more involved due to the double energy integrations in Eq. (351. Details are provided in Appendix [C] 

The second-order correction in Vr to the total current in linear response to the applied bias V is obtained by 
summing and /pj'^ over all scattering states at the Fermi energy. Since the matrix t^iV^'^(£, £)V'^^(£, £) has a 

real trace, the second-order change of the zero-temperature dimensionless conductance read^^ 

g(2) = -4^2 ^r^^ -^12-^21 _ ^,t/ v^iyi^] - 2Re {Tr [tt/V^^V"] }) 

2 



47r ^ Im ■ 



d£ 



(t) £-ei 



Tr[tV'V2'(£F,£)V'l(£,£F) 



(36) 



The first three contributions, where the energy dependence is not explicit, are on-shell terms at e-p, while the last 
contribution involves an energy integration characteristic of second-order perturbation theory. 

Separating in the energy integral the principal-part and delta-function contributions, we can write 



g(2) ^ _ 4^2 (r^^ [^t^ -^12-^21 _ ^,t^, -^21-^12] ^ ^^t^, (^22-^21 _ yaiyll)] |) 

2 roo 



1=1 



(t) £ — £f 



Im {Tr[iV'V2'(£F, £)V'i(£, £f)] } . 



(37) 



As in the case of the first order correction g^^^\ g'-'^^ is 
given by a trace over the incoming right-moving modes, 
and thus independent of the basis. In the transmission 
eigenmodes basis the computation of the trace simpli- 
fies: the terms rH' are relevant outside the regime of 
conductance quantization, but in such the contri- 

bution 5^^^ is negligible (at small Vt) with respect to g'-^\ 
Thus, (7^^^ dominates only on the conductance plateaus 
of a QPC. In this case Eq. (37 1 is considerably simplified 
and becomes 
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(2) 



-47r2 Tr [th V^^V^^] 



N 



= -47r2 y 

Z ^ L -I ? 



(38) 



the sum running over the open eigenchannels of the 
QPC. 



Equations (311 and (37 1 provide all the information 



about the non-invasive effect of an SGM tip applicable 
to situations without and with conductance quantization, 
respectively. While some general conclusions can be di- 
rectly extracted from these expressions, their dependence 
on the unknown matrix elements V''' makes it difficult to 



perform explicit calculations and to address important 
questions like the locality of the response. In the follow- 
ing section we consider the particular case of an abrupt 
QPC since in this system the scattering states can be cal- 
culated explicitly and thus the evaluation of the matrix 
elements becomes possible. 



V. SCATTERING STATES OF AN ABRUPT 
QUANTUM POINT CONTACT 

QPCs constitute a paradigm of quantum transport 
since they lead to the interesting phenomenon of con- 
ductance quantization .'^'^'^ As a consequence, they have 
been intensively studied with the SGM technique .EHii] In 
order to quantify the conductance corrections of Eqs. 



(31 1 and (37 1, we need the scattering states of the un- 



perturbed system. While the conductance quantization 
of a QPC is very robust with respect to the details of 
its geometry,'32Hl3l ^jjg scattering states are highly depen- 
dent on the details of the constriction. In this section 
we focus on the particular case of an abrupt quantum 
point contact (AQPC) with a wide-narrow-wide geome- 
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FIG. 1. (Color online) The wide- narrow- wide geometry repre- 
senting an abrupt QPC. The 2DEG (white area) is delimited 
by gates modeled as hard- wall boundaries (thick black lines). 
The two sketched systems of polar coordinates p, 9 are used 
in Sec. |VI| to express the scattering eigenstates in regions I 
and II. 



try as shown in Fig. [T] In such a setup, the conductance 
was calculated analytically in an approximate way,^ giv- 
ing an understanding of conductance quantization's key 
ingredients. We extend these calculations in order to ob- 
tain the scattering states, and then the matrix elements 

The right-moving scattering states in region I and II of 
an AQPC (see Fig. T]) have the asymptotic form given in 
Eq. (Gal, while within the constriction, i.e. for — L/2 < 
X < Lj2, they are 



(Tie'''"" +7^ae-'^"^) 



with the transverse channel wave-functions in the narrow 
region 



1 



sin [Qniy - w)]. 



(40) 



In analogy with Eq. ^ the transverse energy in the nar- 
row region reads e^*-* — h^Q\/2Mc with the transverse 
momentum Qn ~ 7rn/2w. The longitudinal momentum 
Kn = yjk? — is real for the propagating channels 
and pure imaginary for the closed channels. In the wave- 
function matching at x = ±L/2 the fundamental quanti- 
ties are the overlaps between the a and the n transverse 
channel wave-functions 



where 



VwWl 



(41) 



fn{z) 



cos(z) for n and a odd 
sin(z) for n and a even 



(42) 



By symmetry there is no coupling of channels with dif- 
ferent parity. Notice that the c?„a's are not singular at 
qa — Qni ^nd since 0^ and <!>„ are orthonormal bases 
in their respective y-intervals, we have the completeness 

relation J^a'^-'ia^n'a = Sn7i'- 

There are two important observations at this point: 
(i) the etna's are appreciably different from zero only 
for Qa € AQn = [Qn~i:Qn+i\ and (ii) ka is a smooth 
function of Qa- Based on these observations, Szafer 
and Stone introduced the highly suc cess f ul m ean field 
approximatioEP (MFA), where Eqs. (|4T|)-(|42| are sim- 
plified by taking 



"a I MFA 




if qa e AQn 

if qa i AQn 



(43) 



Such an approximation allows for an analytical determi- 
nation of the conductance of an AQPC, but does not lead 
to the scattering states. Thus it is necessary to have a 
less restrictive approximation, a smooth field approxima- 
tion (SFA) where the qa dependence in iAna of Eq. (41 1 
is kept. 

Next, introduce a generalized momentum-like quantity 
crucial to the wave-function matching problem 



(44) 



The prime over the summation symbol indicates that 
only the a's with the same parity as n are considered. 
When going to the continuum limit in the wide region 
the density of modes q is W/n - i.e. only half the usual 
one. This leads to 



W 

./o 



dqk^'\q)iAn{qyK-{q), (45) 



where fc*-'-* (g) = ^Jk"^ — (p- is the longitudinal wave- vector. 
Taking into account the smoothness condition (ii) and the 
completeness relation of the iJina& gives 



(46) 



Consistently with the above mentioned property (i), we 
take ICn as the average value of the longitudinal wave- 
vector in the interval AQn, that is 



dgfc«(g). 



(47) 



AQ„ 



According to the positioning of k with respect to the 
interval AQn, the generalized wave- vector /C„ may have 
real and/or imaginary parts. 

The diagonal form ( 46 1 greatly simplifies the wave- 



function matching problem, leading to 



Dr. 



2VkJ^{Kn±ICn)i.4 na 

{Kn + Knf e 



-i(fc„±/f„)L/2 



{Kn - ICnf e 



(48) 

(49) 
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The transmission and reflection amplitudes, respectively, 
are 

Kn 

Dn 



tba = 4yfc^e-'[('=^+'=")^/2l J2 jf^nb^^na, (50a) 



Tba 



-2yfc^e-[('=' + '=")^/2l J2 ^^nbAna (50b) 



with 



{Kn +/C„)e 



-iK„L 



+ {Kn - ICn)e 



(51) 



The dimensionless conductance naturally acquires the 
form pO| with 



Tm — 



4|ir„|Re{/C,„} 

I An I 



(52) 



explicitly showing the one-to-one correspondence be- 
tween the transmission eigenmodes m and the constric- 
tion channels n. Therefore, from now on we will make 
the identification n = m. Notice that Eq. ( 52 1 takes into 
account the contribution of the open {Qn < k) and the 
evanescent (Q„ > k) channels of the QPC, since for the 
lowest evanescent mode Re{/C„} 7^ 0. 



The conductance resulting from Eq. ( 52 I is shown in 
Fig. |2] and agrees with that obtained under the more re- 
strictive assumptions of the MFA, providing an extremely 
good approximation for the numerically calculated con- 
ductance of an AQPCl^ll The conductance oscillations 
as a function of k result from quantum interference oc- 
curring in the AQPC.!^ A finite temperature provides 
considerable smoothing and an improvement of the con- 
ductance quantization.'^^ 

The scattering states are completely determined by 
Eqs. (48l-(51), putting us in a position to readily eval- 



uate the matrix elements ( 25 ) once the tip potential is 
specified. In fact, as we saw in Sees. |III| and \rv\ the ex- 
pressions for g^^^ and g^"^^ simplify considerably in the 
transmission eigenmodes basis, which will then be our 
choice in the next section. 



A 



B C 



D 



1 








n 




\ ---7^1 








'1 \ 



3 

kw 



FIG. 2. (Color online) Total transmission (solid) and re- 
flection eigenvalues of the first (dotted) and second (dashed) 
eigenmodes of an AQPC, as a function of the wave-vector k, 
for the lowest 2 channels [obtained from Eq. ( |52[ ) of the SFA] . 
The blue circles indicate the fc-values used in Figs. [3]|6] 



The coefficient cITIn coincides with the matrix element 

1(2), £, a 



[^1(2)] ma '-'f I^Q- ( 19 ' ^^P to an overall m-dependent phase. 

From the form^ ) of the scattering wave-functions and 
the form (50a I of the transmission amplitudes, the gen- 



eral expression for the scattering eigenstates in region II 
is 



Xl e m ^ 



r) 



Dn 



(54) 



The spatial dependence of xYe m through the function 



F„(r) = ^c4„,h (l>b{y) e 



ikh(x-L/2) 



while the prefactor A^n is defined by 



-ikaL/2 



(55) 



(56) 



VI. SCATTERING EIGENSTATES OF AN 
ABRUPT QUANTUM POINT CONTACT 

In the same way as the incoming modes ([2| generate 
the outgoing scattering states (|6| , the transmission eigen- 
modes (19 1 give rise to scattering eigenstates Xi,e.m that 
are eigenfunctions of the current operator ij^^^ ■ Deter- 
mining Qi Jj^^ requires the diagonalization of tH and t'H' , 
which is, in general, a difficult task. The knowledge of 
the scattering states of an AQPC developed in the last 
section will be used to derive approximate expressions for 
the scattering eigenstates. The latter can be written as 
linear combinations of the scattering states 



X,,™(r)=^cS<Jr). 



(53) 



The form 



I kg 

Re{lCm} 



JkaL/2 



satisfies the normalization condition 
leads to 



Ann = VRe{/C„i} S„ 



'I.e. a 



(57) 



1 and 



(58) 



ensuring that xYe m is an eigenvector of I^'^ with an 
associated current eigenvalue 



(0) 



e {A\Krn\)Mmrn\''Re{Krn} 

h \Dm\^ 



(59) 
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Since 

„(™) 



AO) 



Ci J ^ is consistent with the form (52 1 of the transmission 



eigenvalues. The c^™^^ are obtained up to an overall, 
m-dependent phase, as it is usually the case in the deter- 
mination of eigenvectors. Given the simple form ( 58 1 of 
Aran , the wavc-function ( 54 ) of the right-moving scatter- 
ing eigenstate reads 

X"e,m(l-) = , . tm F™(r), (60) 

having defined t„j — (4isrmRe{/C„j})/I?,„. In general 
l^ml = I'm; while for the propagating channels tm = I'm- 
Frai'^) can be calculated by going to the continuum 
limit in the sum over the transverse channels 



|7^j the choice (57 1 for the coefficients 



Fm{^) 



2(-l)™ 
7ri(ii/2 



dqfmiqy) e'" 



»fc("(g)(x-L/2) Qmfmjqw) 



(61) 

Since the function defined in Eq. ( 42 ) is trigonometric 

(62) 



we can split i^m(r) in two integrals and write 
i^„(r)=<;„[F+(r) + (-l)"-iF„(r)] 
with <^jn — I (— i) for odd (even) m and 
1 



has a smooth dependence on q (even for q 
whereas h± is defined as 

h±{q)^k^'Hq){x-L/2)±qy. 



(65) 



Using two-dimensional polar coordinates in region II (I) 
with origin at the center of the right (left) opening 
one can write p — r ^ (L/2)x = p (cos 6*, sin 6*), where 
e [0,7r/2] U [37r/2,27r] for II and G [7r/2,37r/2] 
for I. Expressing the wave-vector in polar coordinates 
{k''^'' {q) , q) — fc(cos?9, sini?) one can write 



(q) — kp [cos ■& cos 9 ± sin d sin 6] 



(66) 



For fcp 3> 1 the functions h± are rapidly varying with 
q (or d), and therefore F^{r) can be computed in the 
stationary phase approximation leading to 



X"e,m(r) = C 

with 

OmiKO) 



^m^ 



-ITT I A 



^Re{A:„,} V TTw ^ 



kQm COS 6 



fm{kw sinO) 



(A: sin 0)2 -g2^ 
Notice that the condition kp ^ 1 means that our results 
hold already at distances from the opening of the order 
of a Fermi wave-length. The asymptotic form of the scat- 
tering eigenstates has the form of a radial outgoing wave 



The function 



TTW 



1/2 



/ dq (g) 6*'*='= (63) e ''/-y/fcp, with an angular modulation 6m (/c, 6*) peaked 

Jo alone the directions 9 for which there is a mode matching 



9miq) 



Ql 



fm{qw) 



(64) 



along the directions 9 for which there is a mode matching 
between the transverse wave- vectors Qm and k\ sin 6*1. 



Using the coefficients c^™^ of Eq. (57 1, a similar 



stationary-phase integration over the transverse channels 
q yields the form of Xi,e,m in region I 



Xl,e,m(l") 



/hp 



-ihp 



Qmik.9) 



(69) 



The refiection amplitude Tm. reads 



2SmRe{/C™} ^ 
Tm = Y) ^ 



(70) 



and is such that 72.^ = Vm^ 1 ^ KnP- The reflection 
and transmission eigenvalues of an AQPC are plotted in 
Fig-i 

Figs. |3] and |4] show, respectively, the scattering eigen- 
states in regions I and II for values of k at the step and 
on a plateau for m = 1 (top panels) and m — 2 (lower 
panels) . The scattering eigenstates show all the expected 
features: (i) a number of lobes dependent on the mode 
index; (ii) a widening conical shape centered on the QPC 
(iii) a decay proportional to 1/ ^J~p. At a conductance step 
the eigenmode is mainly reflected, leading to a highly 



r 



asymmetric wave-function, whereas on a plateau an al- 
most perfect (7^ w 1) transmission yields a rather sym- 
metric pattern (with respect to the y-axis) . As discussed 
in the previous section, the abruptness of the constriction 
degrades the perfect 7^ = 1 transmission'*^ typical of 
a conductance plateau, resulting in the very weak inter- 
ference fringes from reflected waves visible in region I. 

In the limit of a point-like opening, kw <^ 27r, the 
problem becomes analogous to that of two-dimensional 
waves diffracting at a slit, and indeed the wave- function 
reduces to an angle-modulated outgoing Hankel function 
in the asymptotic limit, that is. 



cos{e) 



(71) 



In the opposite classical limit of a wide constriction. 
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?1 




xjw 



xjw 



FIG. 3. (Color online) Color scale representing the squared 
absolute value of the scattering eigenstate wave-functions, 
multiplied by p, in region I [left panels, Eq. ( 69 1] and II [right 



panels, Eq. (67l[, when the energy e is such that the partic- 



ular mode is opening (i.e. at the conductance steps). The 
upper and lower panels, respectively, show p|xi,e,i(r)P and 
p|Xi,£,2(i')P, with a scale given by the extreme values of the 
former. The corresponding value of k considered in the up- 
per (lower) panels is marked as ("C") in Fig. [2] and has 
7?.i ~ 0.2 (7?.2 ~ 0.2). The factor p compensates for the Xj p 
decay of the squared modulus and is introduced in order to 
make the wave-function profile visible at large distances. The 
distances x and y are scaled with the width w of the constric- 
tion, L = 2.5 ui. 




0.5 



xjw 



xjw 



FIG. 4. (Color online) Same as Fig. [s] but for energies e 
such that the particular mode is open (i.e. on the m — 1,2 
conductance plateaus). The k value for the upper (lower) 
panels corresponds to point "B" ("D") from Fig. [5] 



kw ^ 27r, it instead reduces to a straight ray 



Al,£,m 



(r) 



„ifcp 



m 



(72) 



Though x\.e,m(j) was expHcitly computed separately 
in region I and II, in the case of open modes (conductance 
plateaus with 7^ sa 1 and 7?.^,^ « 0) the central symmetry 
of the QPC translates into xi,e,m(r) — Xi e m(~i"); up to 
the above-discussed residual interference fringes linked 
with the non-unitarity of the transmission. Thus, in a 
centrally symmetric constriction without magnetic fields, 
knowledge of the wave-function of an open eigenmode 
impinging from one of the leads (1 or 2) in one of the 
half spaces (a; > or x < 0) is enough to evaluate the 
matrix elements of the perturbation. 

According to Eq. (67 1 the scattering eigenstate Xi,e,m 



has the form of a radial function times a real function 
containing the angular dependence. The corresponding 
current density in region II is a vector field oriented in 
the radial direction p with absolute value 



„-p(o)Ji 



(r) 



ehk 



Al ,e,m 



(r)| 



(73) 



The current density in region I a radial vector field as 
well, whose absolute value reads 



2e 



1 



h Re{]Cm}p 



(-l + |r„p)e2„(fc,0) . (74) 



In a conductance plateau r„ 
proportional to |xL,„(r) 



I = 0, and then jf_'°i;^(r) is 
We stress that such a con- 



nection, like that of Eq. (73), only holds if: (i) the tip 
is at a certain distance from the QPC, pk 1, and 
(ii) the system is perfectly clean. These simple observa- 
tions have a crucial effect on the physical interpretation 
of SGM data, as we discuss in the following Section. 



VII. CONDUCTANCE CORRECTIONS FROM 
A LOCAL TIP IN THE NEIGHBORHOOD OF 
SYMMETRIC AND ASYMMETRIC 
STRUCTURES 

The tip induces a local perturbing potential centered 
on the projection rx of the tip center on the 2DEG. For 
simplicity's sake consider the extreme case of a (5-like 
probe 



VT{r) — VTS{r ~ rx) 



(75) 



as the use of more realistic smoother potential profiles 
would not qualitatively change our conclusions. Accord- 
ing to Eq. (l75| , 



= VT X2,e,m(rT) Xl,e,m(rT) 



(76) 



In general, there is no relcitioiisliip between Xi,£-m 

and 

X2,e,m since the matrices ui and U2 of the polar decom- 



position (17 1, defining the transmission eigenmodes (19 1 
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and the corresponding scattering eigenstates ( 53 1 , are in- 
dependent. Particular cases arise in symmetric setups, 
where the scattering matrix S exhibits further symme- 
tries, beyond the condition S"^ — S stemming from time- 
reversal. Our example of an AQPC is one of them, since 
it has the fourfold (4F) symmetry due to the up-down 
and left-right refection symmetries. The 4F symmetry 
allows to write the scattering matrix ([t]) as 



S = 



5e 




(77) 



where the N x N scattering matrix S'cfo) corresponding 
to even (odd) modes has the structure^ 



c(o) 



'"c(o) ^c(o) 
^c(o) ^c(o) 



with 



c(o) 



re(o) and t^^^ = te(o)- 



(78) 



In the polar de- 



composition of <5'e(o) we see that Ui^o(o) — "^2,0(0) ■ More- 
over, since Eqs. ([2| and (|5| dictate that ip2~J ^{x , y) = 

i-'^YvVai-^^-y) for x > 0, 



X2,, 



.(r)-(-l)'"-W,™(-r). 



(79) 



The first-order conductance correction of a 4F symmetric 
structure perturbed by a delta-function tip potential at 
tt has the form 



N 



(rx) =invTj2 ("1)"'' 



m— 1 

Im{xte,m(-J'T) Xl,e,m(rT)} . (80) 

The dependence of g'^ on the values of the scatter- 
ing eigenstates at two points makes it independent of 
the overall m-dependent phase appearing in the determi- 
nation of Xi,e.m and gives a non-local character to the 
conductance correction. However, since the two points 
involved are symmetric with respect to the center of the 
structure, the corresponding values are indeed related. 

Our example of an abrupt QPC is a particularly simple 
case of a 4F structure. The condition of Eq. (79 1 is read- 
ily checked for an AQPC since X2,s.m can be obtained 
from Xi.e,m by exchanging I by II and r by — r (or 6 by 
e + n). From Eqs. (|67|, and Q, the first-order 



conductance correction is 



kpp I D„ 



(81) 



for an AQPC probed by a delta-function tip when the 
Fermi wave-vector sets the conductance at the m*'' step, 
which is plotted in Fig. |5] 

This expression is valid in both regions I and II, and 
5AQPC,m(-i"T) = 5'AQPC,m(i"T). Such a result is consis- 
tent with the general finding of Sec. Ill that g^^^^ remains 




FIG. 5. (Color online) Conductance corrections at the low- 
est two steps for an AQPC from Eq. (81 1. The 1/p decay is 
compensated as in Figs. [3] and [4] and the color scale is set 
ffAQPC i('"t). The presence of the SGM tip at certain po- 
sitions enhances, rather than lowering, the conductance, in 
agreement with the numerical simulations of Ref . |28) . The k 
value for the upper (lower) panels corresponds to point "A" 
("C") from Fig. [2] 



invariant under the changes Vrix^y) — ?■ VT{±x,±y) 
when the unperturbed system exhibits the correspond- 
ing (up-down and/or left-right) symmetries. 



The first-order conductance correction of an AQPC 
in the m*** step has the same angular dependence as 
the probability density of the corresponding scattering 
eigenstate, oscillates in the radial coordinate as cos2fcF/0, 
and decays spatially as (fcpp)"^. Such a behavior is 
in good agreement with quantum mechanical numerical 
calculations^^. Even if 3AQPCm('"T) appears as a local 
correction, it is not directly related with local observables 
such as the probability density or the current density. 
Moreover, we stress that the local character of Saqpc m 
stems from the symmetry of the unperturbed structure. 
For a non-symmetric structure the first-order conduc- 
tance correction is in general not local, and unrelated 
to the probability density or current density. 

On a conductance plateau the dominant correction is 
f;(^^ and, according to Eq. (38 1, for a (5-function tip at tt 



N 



(82) 



In a 4F symmetric QPC the condition (79), together 



with the relationship between the probability density and 



the radial component of the current density ( 73 ) , yields 
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FIG. 6. (Color online) Same as for Fig. [sl but on the lowest 

J V — f 2) 

two plateaus, Eq. (84 1. The color scale is set by ffAQPc(''T) 



on the first plateau. The k value for the upper (lower) panels 
corresponds to point "_B" ("D") from Fig. [5] The minima of 
the conductance correction correspond to the maxima of the 
squared absolute scattering eigenmode wave- function. 



in both regions I and II, 



f4F-QPC 



Me 



N 



(83) 



th 



where fcp is such that the conductance is on the N' 
plateau. The second-order correction in this case is then 
local, and proportional to the square of the total current 
density. 

For the case of an AQPC the conductance correction 
on the N^^ plateau takes the simple form 



5aqpc(i"t) 



2M, 



N 



TTwh? kpp 



Re{/C„ 



} 



(84) 



and is plotted in Fig. |6] Notice that our expression de- 
rives from Eq. ( 38 1 . Thus we are approximating the full 
expression ([37]) for g'^^ by taking — 1, and therefore 
neglecting the small reflection visible in the eigenmode 
structure of Fig. |4] As in Eq. (811, the above expres- 
sion is symmetric, and valid in both regions I and II. All 
open channels contribute and the overall radial decay is 
as (fcpp)"^, consistently with the quantum mechanical 
numerical calculations .f^Sl It is important to remark that 



the simple connection ( 83 ) between the second-order con- 



ductance correction and the current density is not univer- 
sal, but only holds under the very restrictive conditions 
discussed above. Even an ideal probe would not be capa- 
ble of imaging the current density in a system that does 
not show conductance quantization or if the latter arises 
in a non-symmetric structure. 



While we expect that realistic QPCs can be approxi- 
mately symmetric,^ it is clear that once disorder is intro- 
duced in the 2DEG neighboring the QPC, the symmetry 
of the unperturbed structure would be completely lost on 
scales larger than the elastic mean-free-path. 

Once the 4F symmetry is broken the connection be- 
tween the two scattering eigenstates xi,e.m and X2,e,m is 
lost, since the matrices ui and U2 of the polar decomposi- 



tion (17) are independent, and therefore the conductance 



correction (82) depends on the probability densities as- 
sociated with two different states. Moreover, in the case 
where there is scattering between the QPC and the tip, 
the separability of the scattering eigenstate in the radial 
and angular coordinates is lost, leading to a current den- 
sity which is no longer directed along the radial direction 
and proportional to the probability density [as in Eqs. 



(731 and (74l] 



From the above considerations we conclude that a 
quantitative analysis establishing the regimes for which 
the proportionality between the conductance correction, 
due to the presence of a weakly invasive tip, and the 
square of the current density holds, crucially depends on 
the symmetries of the unperturbed structure and on the 
potential landscape generated by the neighboring impu- 
rities of the 2DEG. The relation between the measured 
conductance change and the current density is not univer- 
sal, but depends on the particular SGM setup. Though a 
more detailed study of the scattering eigenstate profiles 
is beyond the scope of our present work, it is a relevant 
issue for a conclusive, clear-cut interpretation of SGM 
experiments performed in generic geometries and under 
different regimes. 



VIII. CONCLUSIONS 

In this work we have presented the theory of scan- 
ning gate microscopy applicable when a weak probe is 
scanned above a 2DEG in the neighborhood of a nanos- 
tructure. We have discussed and extended the pertur- 
bative approach put forward in Ref. [28. and applied it 
to the paradigmatic case in which the nanostructure is 
a QPC. In the case of an AQPC an approximate solu- 
tion of the scattering eigenstates is derived and used to 
understand the effect of a weakly invasive SGM tip. 

The first-order conductance correction in the tip po- 
tential involves the unperturbed scattering amplitudes 
and scattering states impinging on the structure from 
opposite sides. In the case of a QPC this first order is 
the dominant correction only in the conductance steps of 
partially open eigenmodes, but it is suppressed on con- 
ductance plateaus. This correction is local in space, os- 
cillates on the scale 2kp, and decays as the inverse of the 
distance between the QPC and the tip, but it is not di- 
rectly related with the probability density of the scatter- 
ing eigenstates, nor with their current density. When the 
Fermi energy is set on a conductance plateau, the second 
order correction in the tip potential becomes dominant. 
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In the case of a symmetric structure it becomes a local 
quantity, given by the square of the local current density 
arising from all the open eigenmodes. 

The dependence of the conductance correction on tip 
voltage has been observed to be linear in the case of a 
non-quantized conductanc^^Sl^ jn agreement with the re- 



sult (|31[) for g'^^\ The quadratic dependence of g^^^ on 



(2) 



the tip-induced potential Vt is an important prediction 
that should be observable on the conductance plateaus. 

Within the framework of linear response theory in the 
applied bias voltage V, the detailed local information 
encoded in the current density and the electric field is 
irrelevant .SSI Only integrated quantities, the total cur- 
rent carried by the scattering states at the Fermi energy 
and the applied voltage, appear in the determination of 
the conductance. The spatial resolution provided by the 
SGM setup yields conductance corrections that could be 
interpreted as a measure of local properties. However, 
even in the least invasive case, it is only under very spe- 
cific conditions of symmetry and unitary propagation of 
the eigenmodes that these corrections could be shown to 
be local, and related to the current density. 

The abrupt quantum point contact is a good m odel 
for the constrictions used in some of the experiments.l^ESl 
Constrictions with more gradual openings exhibit colli- 
mation phenomena, where the momentum distribution 
of the electrons leaving the junction is weighted towards 
large longitudinal momentum.!^ In such cases, the cur- 
rent densities are more forward focused than those pre- 
sented in Figs. [3] and |4j The corresponding experimen- 
tally observed conductance corrections also exhibit some 
forward focusing .'33 Even though an analytic determina- 
tion of the scattering eigenstates is not possible for such 
geometries, the qualitative conclusions of our calculations 
remain valid, and the symmetry requirements necessary 
for establishing a connection between the conductance 
correction and the local current density also hold in this 
case. 

There are a number of ingredients, beyond the scope 
of the present paper, that might be of quantitative rele- 
vance in certain parameter regimes. Among them, non- 
perturbative effects in the tip potential, non-linear effects 
in the bias voltage, as well as interaction, temperature, 
disorder, and magnetic field effects. 

The tip potentials currently used are strong enough to 
create a divot of depletion in the 2DEG,'^ and multiple 
scattering between the nanostructure and the tip. This 
is not a serious restriction for the perturbative approach 
since its applicability relies on the smallness of the ma- 



trix elements ( 25 1 , which become progressively smaller 



for large distances between the tip and the nanostruc- 
ture. Moreover, our conclusions on the necessary con- 
ditions for establishing a connection between the correc- 
tions and the unperturbed current density, in the least 
invasive case, are expected to be more stringent in the 
non-perturbative regime. The Green function approach 
sketched in Appendix [D] is cumbersome due to the pro- 
liferation of energy integrals, but it is amenable to be 



extended to higher-order in the perturbation scheme. 

While most of the SGM experiments have measured 
linear conductances, some recent experiments investi- 
gated the non- linear regime.'^ The non-linear theory of 
SGM is considerably more complicated than its linear 
counterpart described in this work since electron-electron 
interactions have to be incorporated at least in a self- 
consistent form in order to respect charge conservation. 
The spatial symmetry breaking induced by the perturb- 
ing tip might give rise, under certain conditions, to a 
measurable asymmetry of the I-V characteristic. B21 

The electron-electron interactions are crucial in the 
non-linear regime, but they are also expected to be im- 
portant in the linear one, in particular in the case of 
partially open modes.^^S The Friedel density oscillations 
induced by the tip are expected to modify the electro- 
static potential in the constriction. This effect, different 
from the cross talk between the tip and the gates defin- 
ing the structure,'^ has been treated in one-dimensional 
lattice models in Ref . HHl 

Temperature effects are important to suppress quan- 
tum interference within the constriction and to improve 
the conductance quantization."^^ When the smoothing 
provided by a finite temperature is large enough, the step 
and plateau behaviors get mixed, and non-monotonous 
dependencies of the fringes appearing in the conductance 
corrections on temperature can be obtained.'^Sl 

Disorder is always present in the 2DEG neighboring 
the nanostructure. The sample quality translates in the 
disorder strength of the 2DEG measured in terms of the 
mean free path, and the dominant scattering characteris- 
tics of short- or long-range disorder.^ ^ When long-range 
potentials leading to small-angle scattering are dominant, 
the SGM sweeps reveal a branching structure that has 
been interpreted in terms of caustics related with the 
classically underlying electron motion. The presence of 
disorder in the 2DEG breaks any possible spatial symme- 
try of the nanostructure under study. It remains to be 
determined beyond which distance the symmetry break- 
ing becomes relevant. 

The SGM setup has been used to study the focus- 
ing effect of a perpendicular mag netic field.l^ The im- 
portance of interference between paths traveling between 
two nearby QPCs, with and without defiection from the 
tip has been demonstrated. The theory presented in the 
present work can be easily extended to include the time- 
reversal symmetry-breaking effect of a magnetic field. 



Appendix A: Energy integrations for the current 
corrections 



In this appendix we present the energy integrations 
appearing in Eqs. (26 1 and (CI I, which are characteristic 
of scattering theory in wave-guides .I2il3l! The general form 
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Im Ikr, 



-hi - in 



ka + ir] 



Ro kr, 



it 



FIG. 7. (Color online) Contours in the complex ka plane for 
the integration of Eq. ( A2 1. The poles shown are for an open 
channel a. For sx > (< 0) the contour ri(r2) has to be 
chosen. This way the integral on the quarter circle (dashed 
lines) goes to zero for fc^ — > oo, whereas the one along the 
imaginary axis (dotted lines) becomes Laplace-like and thus 
0{l/\x\). 



of these integrals can be expressed as 



(t) 



N 



(Al) 



where fcg = y 2Afc(e — £a*')/^^- The infinitesimal pos- 
itive (negative) part given to applies to terms aris- 
ing from outgoing (incoming) modes, while s = ±1 ac- 
cording to whether the corresponding mode is a right or 
left-moving one. In Sees. |III| and |IV| we evaluate cur- 
rents in the right lead, we thus are interested in the case 
i/2 ^ a; < oo, and the sign of sx is that of s. In Ap- 
pendix [D] we have the opposite case since a;,x < 0. The 
function ^a'a(£^j£) is assumed to have a smooth depen- 
dence on e and e. 

Changing the integration variable from e to fcg we have 



-2E 



\ 2Afo 



(t) 



[ka - {ka + ir])] [ka + (ka + iri)] 



(A2) 

where fcg — \J 2Mc{e — e"^^)/'h? is real for an open chan- 
nel a and imaginary for a closed one. For an open channel 
(a < N) the poles in the complex fcg plane are indicated 
in Fig. [7] The integral along the positive real axis can be 
obtained from the contour integration along the paths Fi 
and F2. The choice of the appropriate contour depends 
on the sign of s in order to ensure that (i) the contribution 
from the quarter circle (dashed lines) vanishes exponen- 
tially in the limit |fca| — > 00; and (ii) the contribution 
along the imaginary axis (dotted lines) is asymptotically 
small, 0(l/|x|) in the limit x — >■ +00. For sx > 0, Fi has 



to be used, and the integral is given by the residue at the 
pole ka = ka + irj. For sx < 0, F2 has to be used, leading 
to a vanishing integral. For a closed channel (a > N) the 
poles are close to the imaginary axis and both contours 
lead to a vanishingly small contribution. We then write 



N 
5=1 



(A3) 



The first term of Eq. (27b) has the form above described. 



and the corresponding e integration leads to a vanishing 
contribution. The contribution from (27a I and the sec- 



ond term of Eq. (27b) are slightly different, since they 



can be written in the form 



de 



N N 



(t) 



5^5]el(e,e,&)e-^?% (A4) 



a=l 6=1 



with N = min{A^, A^}. Now, it is convenient to pass 

from the variable fcg to fcf, = 
and write 



N 



3=1 6=1 



iaa (S +4*^£.& 

[^6 - (^6 + iri)] [h + ih + ifl)] 

(A5) 



where kb = \J 2Me{e — e'^^) /h? is always real since h < 
N. The lower integration limit fc™™ is equal to if a < 

b, and 2Mc (e^ "* — e^'"* ) /K^ if a > 6. In the first case 

(only occurring when a < N) the pole at k^ — k^ + if] 
is inside the contour Fi of the complex fcf, plane. In the 
second case the contour Fi has its vertical component 
with Re{fcb} = fc™'" and the pole is within Fi only if 

Eg*-* < e, that is, if a < A^. The resulting non-vanishing 



contributions to Eq. (A4| are 



N N 



(A6) 



a=l 6=1 

The application of Eq. (|A6|to the e integration of 



Eq. (26) readily leads to Eqs7l28|-(29) 




Appendix B: First-order conductance corrections 
from functional derivatives 



In this appendix we show that the first-order correction 
(31) contains as a particular case, the conductance 
301 



sensitivity (obtained in Ref. '30| to a change of the scat- 
tering potential U{x) defining a one-dimensional struc- 
ture. One-dimensional systems are particularly simple 
within our approach because all the scattering states are 
scattering eigenstates. 

Writing the dimensionless conductance as g = l^tPj 
where tt is the total transmission amplitude taking into 
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account the unperturbed scatterer and the perturbation, 
the variation of the conductance when the scattering po- 
tential is changed from U{x) to U{x) + 6U{x) is 



Sg = —Air / dxr]{x)SU{x) 



(Bl) 



The function r/ is, up to a niultipHcative constant, the 
functional derivative of the conductance with respect to 
the electrostatic potential, and can be written as 



with all quantities taken at the Fermi energy. The trans- 
mission and reflection amplitudes of the unperturbed 
delta-barrier potential are t = t' = 1/(1 + iu) and 
r = r' = — m/(l + iu), respectively. Therefore, 

^-AttuT I d2;Re{*^°'*(x)*f'(a;)}(5C/(x). (B5) 



The expressions (|6| of the lead states yield 



r]{x) = - 



1 

47r 



5U 



St: 



SU{x) SU{x) 



tA . (B2) Re|*^"^*(a;)*^°^(a;)| = ^ [cos(2fcFa;) -|-usin(2fcF|a; 



The case where the unperturbed scattering potential is 
a delta-function U{x) = u6{x) allows for an explicit eval- 
uation of the functional derivatives. The unperturbed 
transmission probability for a particle at the Fermi en- 
ergy is given by T = 1/(1 + with u = u/{hvp), where 
vp is the Fermi velocity. The sensitivity ri(x) is given 



r]{x) 



hvp 



[cos(2A;Fa;) -t- 'usin(2fcF|x| 



(B3) 



We now approach the problem from the perturbative 



expression (31 1 that for a one-dimensional system under 
a change SUjx) of the scattering potential yields the first- 
order conductance correction 



47rlmir*i' / dx¥^^*ix) SU{x)¥°\x) 



(B4) 



By plugging Eq. (B6l in Eq. (B5l, one can see that the 



first-order conductance correction obtained from our per- 
turbative approach can be written in the form Eq. (Bl I, 
with a sensitivity ri{x) that coincides with the result (B3 1 
from the functional derivative approach of Ref. |30| . 



Appendix C: Second-order correction to the current 
carried by a scattering state 

In this appendix we provide the details of the calcula- 
tion of the second-order correction, in the tip potential 
Vt, of the current carried by the scattering state l,e,a. 

with 7 = a, l3, of the 



The two contributions 
second-order current density corrections ( 35 1 , when inte 



grated over a cross section Sx at the right of the scatterer, 
lead to the corresponding changes of the current 



eh 


Ll=l 


r°° de 


r°° ds 


M 


X(t) e+-e 


J,w e+-e 


eh 


±Re< 

1,1=1 


' de 


df 


M 




X(t) e+ - f 



N N 



^^Zi|,(e,e) [V^t^{e,E)[VT^]'l{i,e) 



a=l a=l 
N N 

EEi 

a=l a=l 



(Cla) 
(Clb) 
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Here, is defined in (27) and 



N 



6=1 



Klie, e-)^c'J2^^^ tlJ,, exp [*(fc+ - k-)x] , 



6=1 



N 



b=l 



Xll{e,e) = a I -Sa, exp [-z(fcr - k^jx] + 



f~ "ha 

kb 

— r 



ha exp - fcj )x\ ^ , 

exp [i{k+ + k+)x] 



N 



r'-aa exp + K )^] + 



'bVb-a exp \i{k+ - k^ )x] ) , 



6=1 



(C2a) 
(C2b) 
(C2c) 



(C2d) 



where the barred longitudinal wave-vectors and scat- 
tering amplitudes are taken at the total energy e, and 
N = min{7V, N}. 

For the current correction I^l°^^{x), only one of the 
energy integrations can be done as in Appendix [X] For 
I^l^a{^) t^e e' integration can be done as in Appendix [a| 



I 

while the e integration is analogous, up to the fact that 
the poles are in the second and fourth quadrants of the 
complex plane for fcg and fcf,. We then have 



r(2)7 
'l,e,a 



r(2)7 

-1,£ 



(C3) 



with 7 = a, /3 and the N x N matrices 



1>£ 



r(2)/3 



2e 



El-' 



;=i 



de 



(t^tyi' (e, e) (£- e) + t V'V^' (e, e) V'^ (e, e)) 



27rfi 



Re { V" (e, £)tnv" (e, e) + V^' (e, e)r' (e, e) + 2V" (e, £)tV' (e, e) } 



(C4a) 
(C4b) 



The corrections (C3), summed over a, lead to the second 



order correction given in Eq. ( 36 ) of Sec 



IV 



Appendix D: Green function approach to 
conductance corrections 



In this appendix we provide an alternative derivation 



of our main results Eqs. (31 1 and (36) based on the per- 
turbative expansion of the total Green function. Dyson's 
equation 

g(r,f,£) -g(")(r,r,e)-f J dr g^''\r,r,e)VT{r)g{v,v,e) 

(Dl) 

yields an equivalent approach, with respect to the 



Lippmann-Schwinger equation (21), to the perturbative 



treatment of the scattering problem. 

The Landauer-Biittiker equation (|16| for the unper- 



transmission amplitudes and the Green function, can also 
be used to relate the total conductance g with the total 
transmission amplitude t^, and the total scattering am- 
plitudes in terms of the total Green function Q. Going to 
second order in Vt in the expansion of the total reflection 
amplitude rt ~ r + SA^^ + (Jr*-^-* we have 



(Jg^i) = -2Re{Tr [rUr^^)] } , 



with 



-2Re|Tr 
-Re |Tr 



}■ 



(D2a) 
(D2b) 

(D3a) 
(D3b) 



where the energy arguments are taken at ep. Using the 



turbed problem and the relationship (^9| between the first-order approximation of Eq. (Dll one gets 
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frW ^ in{vavi,y/^ exp [-z(fc+T + k+x)] dy dy / dr 0^(y) g^°Hv,v,e) FtI?) g^''Hr,r,e) Mv) ■ (D4) 



From the spectral decomposition (22 1 of C/'^*'^ in the basis of scattering states the definition (25 1 of the matrix 

elements of the perturbation, and integrating over the transverse coordinates y and y we have 

df f°° de 



<5.« = ^^^^^ expH(fc+x + fc+^)] 



2tt 



1 



TV 



^ {((^aaexp 

a,a— 1 

(^(5aa exp 
taa exp 



^aa exp 



(Sba exp [ifc^ x] + rfca exp [-ifc^x] ) [Vt]^= (e, e) + 



''aa exp 



77* 

iaa exp 



ik~x 



^'te exp [-ifc^x] [yT]aa(^'^) + 
(Jba exp [ifcf^a;] + na exp [-ifc^^a;] ) [VrJal (e, f) + 
fl,exp[-*fc+x] [yT]i(e-f)} . 



(D5) 



Performing the e and e integrations according to the with the matrices defined as in Sec. |III| from Eq. ( |25 1 
prescription of Appendix \K\ for the case x,x < 0, the As noticed in Sec. |III| the first term gives a vanis 
previous expression reduces to 



contribution, and therefore we recover Eq. (31) 



amg 



(D6) 



According to Eq. ( |D2a ) 



g{i) ^ {Tr [rV V" + rh' V^^] } , (D7) is 



The second-order variation of the reflection amplitude 



(pl{y) Cy(°)(r,r,e) V^{f) g('Hr,f,e) V^{f) C/(°)(f,f,e) 0,(y) 



(D8) 



The structure of (|D8| is rather similar to that of (|D4|) for the correction (5r^"'^-'. The main difference is that the 



spectral decomposition of the unperturbed Green function now leads to three energy integrations, where two of them 
can be done as before, yielding 



2 „oo 1- iv 



1^1 



a,a—l 



Then, 
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{2)a 



-47r ^ I TrTr 



1=1 



Im <! r Tr 



rh' V2'(eF,£)V'i(£,eF) 



where the non-specified energy arguments are understood to be taken at ep- Since 

g(2)/3 = _47r2Tr [r^r V"V" + t'h' V^^V^^ 2Re {rh' V^^V"}] , 
the second-order conductance correction can be written as 

^ _4^2 (r^^. ^^,t^, -^21-^12 _ -^12-^21] _ 2Re {Tr [rt^'V^iy"] }) 

de 



1=1 



47ryim<^ / — ^ Tr[rt<'v2'(eF,e)V'i(e,eF) 



(D9) 



(DIO) 



(Dll) 



(D12) 



Using the unitarity condition for 5 we see that Eq. (D12) is equivalent to Eq. (36 I. The symmetry between 



19 



the two ways of writing g^"^^ arises from the different 
choice of the cross sections (a;, a; < in this appendix, 
and X > in Sec. IIVI). 



Appendix E: Exact results in one dimensional 
models 



dimensional system is 



gid = 



1 



1 - 2u Im {r'e*"} + {l + \r'\^ + 2Re {r'e*"}) ' 

(E2) 

We note a — 2kx'T and v — v^/hv. The first and second 
order conductance corrections should then be 



In this appendix we treat the exactly solvable case of 
a one-dimensional scatterer perturbed by a local tip, and 
we compare the resulting conductance corrections with 
those predicted by Eqs. (311 and (37). We consider a 



scatterer described by a 2 x 2 5 matrix, with the form 
([?]) , under a perturbation 



V't{x) — vtS{x — xt) 



where xt is the distance between the scatterer and the 
tip. The dimensionless conductance of the combined one- 



(E3a) 



(2) ~2 

9id = -V 



\t\^(l 



„'|2 



2Re{r'e'"}) -4 [Imjr'e^"}]^ 
(E3b) 



In order to make the connection with Eqs. (31 1 and 



(El) (37 1 the matrix elements (25) involving the unperturbed 



scattering states and the tip potential (Ell need to be 



evaluated. This is easily done in the one-dimensional 
case that we are treating 



V"(e-£) 



^ Y = i*i exp [i{k - k)xT] , 

— a/ = t exp [ifcxx] jexp rifcxxl + r'* exp T— ifcxxl | , 
27r V fc 



27r 



jexp [ifexx] + r'* exp [— ifexx] } {exp [—ikxT] + r' exp [ifca;T]} 



(E4a) 
(E4b) 
(E4c) 



The first order correction (31 ) only involves matrix elements evaluated at the Fermi energy ep, and can be written as 

= -47r Im{r*t' V^^} = 2v \t\^ Im{r'e''=^^ (g^fcxT _,. /*g-.fc:.T) } ^ (E5) 



which is equivalent to Eq. (E3a). Using the form (E4| of the matrix elements, the last contribution of g'-^^ in Eq. (36) 
can be written as 



2d 



Imi (|rf +r'e2'^-"^) 



de 



£ — e 



1 + Re 



2^2 |t|2 (|r'|2(i + Re {/e^"}) + Re {r'e'"} -f Re {(/e^")^}) 



where all terms are evaluated at e in the final expression. When combined with the on-shell terms of Eq. (36) the 
second-order correction of the one-dimensional model becomes 



9 



(2) 



|tp(l 



„'|2 



2Re{r'e'"} -t-2(Re{r'e'"})2^ 



(E6) 



which is easy to verify that coincides with Eq. ( E3b I . 
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